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Different formalisms for unintegrated parton densities are discussed, and some 
results and applications are presented. 

1. Introduction 

Calculations based on coUinear factorization work very well in many appli- 
cations, including e.g.: 

- DIS at large Q^. 

- High p± jets in pp collisions. 

- Inclusive observables. 

In these cases DGLAP evolution works, and fcj^-ordered chains up to a hard 
subcollision or a highly virtual photon dominate the parton evolution. The 
unintegrated parton density, J^{x, fc^), then satisfies the relations 



■ (1) 



This formalism has, however, problems for observables which are sensi- 
tive to the k± of the "active" quark. Some examples are: 

- Transverse momentum unbalance in 2-jet events. 

- Heavy quark production. 

- Forward jets at small XBj- 

In many of these cases calculations, which allow for one extra gluon, e.g. 
LO pQCD + parton showers or NLO DGLAP calculations, are able to give 
a good description of the data, but also these calculations have problems for 
observables which involve a large rapidity separation. In the following I will 
discuss different formalisms for /c^-factorization, non-fc^-ordered evolution 
and some applications. 
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2. Non-fej_-ordered Evolution and fej_-Factorization 

At small X and limited non-fcx-ordcred chains give important contri- 
butions. In a formalism based on fc_L -factorization the non-integrated pdfs 
J^{x, k'j^,Q^) may depend on 2 scales, the k± of the active parton and the 
virtuality of the photon or the hard scattering. The second scale is then 
related to a limiting angle for the emissions, as discussed further below. It 
is also important to remember that the gluon distribution is not an observ- 
able, and depends on the theoretical formalism. 

2.1. BFKL 

In the BFKL evolution^, acciuate to leading logl/.T, J-{x,k'^) depends 
only on a single scale k'j_, and satisfies an integral equation with a kernel 
^(k±,ki): 

J^{x,kl)^To + ^ [- [ d?k'^K{\iMT{-,kl). (2) 

TT J Z J Z 

We note that the dominant leading log behaviour originates from the 1/z 
pole in the splitting function. This leading order result has problems be- 
cause the NLO corrections are very large. There are e.g. large; effects from 
energy conservation, as demonstrated e.g. in MC studies by J. Andersen et 
al? , which imply that analytic calculations often are unreliable. 

The kernel K in Eq. (2) describes the emission of a quasi-real gluon 
with transverse momentum from a virtual link with momentum k'^, 
which after the emission gets momentum kj^ = k^ — qj_. The kernel has 
the property that small values of q_L are suppressed. Such emissions are 
compensated by virtual corrections, and in the BFKL formalism this is 
taken into account by treating the links as Reggeized gluons. 

We can compare this situation with e^e~-annihilation. Here the total 
cross section is determined by the lowest order diagram (a^). The lowest 
order contribution to 3-jet events is ©(as), and to this order there are 
negative contributions to (T2jet5 such that cTtot is approximately unchanged. 
In this case the gluon emissions can be treated with Sudakov form factors. 

For a link in the BFKL chain the 0{as) corrections give a compensa- 
tion for emissions for which < kj_, k'^. These soft emissions give no 
contribution to the inclusive cross section, i.e. to F2. They must, how- 
ever, be added for exclusive final states, with appropriate Sudakov form 
factors. The net result of this is that downward steps in k± are suppressed 
by a factor k\/k^. In the relevant variable hik\, this corresponds to an 
exponential suppression allowing downward steps within ~ 1 unit in \n.k\^ . 
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2.2. The CCFM model 

An evolution equation which interpolates between BFKL and DGLAP was 
formulated by Catani, Ciafaloni, Fiorani, and Marchesini''. This formalism, 
the CCFM model, is based on a different separation between initial and final 
state radiation (denoted ISR and FSR respectively). Some soft emissions 
are included in the ISR, but this increase is compensated by "non-eikonal" 
form factors. The ISR ladder satisfies the following constraint: 

In the ISR emissions the colour order agrees with the order in energy 
(or order in lightcone momentum p+ = Po + Pl)- (Thus all emissions 
which in colour order are followed by a more energetic one are treated as 
FSR.) With this constraint colour coherence implies that this ordering also 

coincides with the ordering in angle, or rapidity. 

This angular ordering implies that the emission angle for the last emis- 
sion must be known, when a new step is taken in the evolution. This angle 
therefore appears as a limiting angle in the non-integrated distribution func- 
tion used in the evolution. Thus the distribution J-'ix, fc^, q) depends on two 
scales, k"^ and q, where the limiting angle is specified by yumit = ^n{xMp/q) 
(in the proton rest frame). 

For exclusive final states, the final state radiation has to be added in 
appropriate kinematical regions. As the initial state radiation is ordered in 
P-^ but not in p_ , also the final state radiation is unsymmetric with respect 
to the initial proton and photon directions. 

We can compare with the BFKL and DGLAP formalisms, where the 
parton distributions depend on a single scale. The initial state radiation is 
strongly ordered either in q+ (for BFKL), or in q± (for DGLAP). In both 
cases this ordering also implies a strong ordering in y. CCFM interpolates 
between the two regions at the cost of a more complicated formalism. 

That more emissions are included in the ISR implies that the average 
step is shorter. Therefore the angular ordering constraint becomes more 
important and implies a strong dependence upon the scale q. This feature 
also implies that it has not been easy to implement the CCFM model in an 
event generator, but such a program, Cascade by Hannes Jung^, is now 
available. The CCFM model does not include quark links in the parton 
chains. The original model, and also the first version of the Cascade 
program, referred to as JS, also included only the singular term oc 1/z 
in the splitting function. To include the non-singular terms is not straight 
forward, but one possible solution is implemented in a new fit, called J2003^. 
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In this fit set 1 includes only the singular terms, while in set 2 also the non- 
singular terms in the splitting function are included. 



2.3. The Linked Dipole Chain Model 

The Linked Dipole Chain Model, LDC'^, is a reformulation and generali- 
sation of the CCFM model. It is based on a different separation between 
initial and final state radiation. The ISR is ordered in both and q-, 
and satisfies the constraint q±i > min(fc_L,, Softer emissions are 

treated as final state emissions. In this respect the LDC formalism is more 
similar to BFKL. Wo note that the ordering in q^ and (j_ also implies an 
ordering in angle. An important property is also that the parton chain is 
fully symmetric with respect to the two ends of the chain. 

The fact that fewer emissions arc treated as ISR implies that a sin- 
gle chain in LDC corresponds to the collective contributions from several 
possible chains in CCFM. Now it turns out that summing over all possible 
emissions in CCFM, the non-eikonal form factors exactly cancel. The result 
is a simple evolution equation in terms of a single scale unintegrated density 
function jF(a;, /c^). In the MC implementation it is, however, also possible 
to add an angular cut and thus obtain results for a two-scale distribution. 

We note that to leading order the LDC and CCFM formalisms give the 
same result for the integrated structure functions. The parton chains and 
the unintegrated distributions differ, however, and only after addition of 
final state emissions in the different relevant kinematic regions do the two 
formalisms also give the same result for exclusive final parton states (to 
leading order). 

Thus the LDC formalism results in a much simplified evolution equation. 
Other merits of the formalism include: 

- It contains the same chains as in DGLAP for large, which makes 
it easier to interprete the differences between large and small Q^. 

- There is a natural generalization to include subleading terms, e.g. 
quark links, non-singular terms in the splitting functions, and a running 
as- 

- It is suitable for implementation in a MC, and the event generator 
LDCMC is produced by Lonnblad and Kharraziha^. 

- The MC gives very good fits to experimental F2 data, and it also agrees 
well with MRST and CTEQ results for the integrated gluon distribution^. 
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2.4. Other formalisms for unintegrated parton densities 

A different formalism, which also interpolates smoothly between BFKL and 
DGLAP, was formulated by Kwiecinski, Martin, and Stasto (KMS)i°. This 
formulation is based on a single scale evolution equation. The contribu- 
tions from > are neglected, and thus the gluon distribution satisfies 
Eq. (1). This formalism has been further developed by Kimber, Martin 
and Ryskin (KMR)^^. In their approach the single scale KMS evolution 
is used, but an angular constraint is applied for the last step. The result 
is therefore a density distribution, which depends on two scales. Since the 
underlying KMS evolution does not include the soft initial state radiation 
in the CCFM model, it also takes fewer and larger steps, which implies that 
the dependence on the limiting angle is significantly smaller than for the 
CCFM formalism. 



3. Comparison between results from different formalisms 

We here want to compare the parton densities obtained in different for- 
malisms, and also study the effects of the non-singular terms in the splitting 
functions and of quark links in the evolution chains. Some results are shown 
in Fig. 1^ Among the LDC fits the result denoted standard contains both 
quark links and non-singular splitting terms, gluonic contains only gluon 
links, and leading only singular terms. (In the fit gluonic-2, refered to in 
the next section, the power of (1 — x) in the input gluon density is changed 
to 7 instead of its value 4 in gluonic.) In all cases the input distributions 
for fc^ — Qq are adjusted to give good fits to experimental F2 data. The 
CCFM result .12003 set 1 contains only gluons and only the 1/z pole in the 
splitting function. We see that the effect of quarks and non-singular terms 
become larger for larger k±. The sensitivity to the scale q is illustrated in 
Fig. 2a, which shows the gluon density J^{x,k^,q) as function of q/k± for 
fixed k±. As discussed above, in the CCFM formalism the density is very 
sensitive to the q scale, and varies strongly for q/k± between 1 and 2. 

It is, however, interesting to note that observable cross sections differ 
much less than the parton densities, as seen in the next section. The reason 
is that the distributions differ in particular for q < k±, while in a hard 
subcollision the dominant contributions are obtained for fc^ « g^/4 ^. In 
Fig. 2b we see that the results are indeed not so different for q = 2k±. In 
this figure also the single scale result from KMS and the derivative of the 
GRV result (denoted dCRV) are included. 
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Figure 1. Comparison of different sets of unintegrated gluon densities at scale q = 10 
GeV. For the notations, see the main text. 
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Figure 2. The LDC gluonic unintegrated gluon density compared to the results of JS 
and KMR. Left: As function of q/kj_ for fixed x and kj_. Right: As function of x for 
fc^ = 10 GeV'^ and q = 2k± . Also shown are here the single scale results from KMS and 
the derivative of the GRV fit. 



4. Applications 

4.1. Heavy quark production 

Applications of the A;^ -factorization formalism to b quark production at 
the Tevatron are shown in Fig. 3. We see that for both the LDC and 
the CCFM models the fits with only the leading term reproduce the data 
best, while the other fits, which should be expected to be more accurate, 
although being significantly better than the NLO QCD result, do not give 
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equally good fits. Wc should note, however, that b production is also well 
reproduced by collinear factorization plus parton showers, as implemented 
in the Pythia event generator ^'^ 




modeli". Right: DO data^^ compared with CCFM results^^. Here also a NLO QCD 
result is included. For notations, see the main text. 



4.2. Forward jets 

Results for forward jet production at HERA are shown in Fig. 4. Here we 
have a comparatively large separation in rapidity, and we see that LO and 
NLO dijet calculations arc far below the data. The CCFM model gives a 
much better description, but also here we see that the best fit is obtained 
including only the singular terms in the splitting function. The reason 
for this is still not understood; is there some dynamical mechanism which 
somehow compensates the effect of the non-singular terms? 

4.3. Minimum bias and underlying events in pp collisions 

In hadron-hadron collisions collinear factorization works well for calcula- 
tions of high-p_L jets. However, in this formalism the minijet cross section 

diverges with a-jet ^ I/pIj which implies that also the total E± diverges. 
This implies the need for a soft cutoff, and in Pythia fits to experimental 
data give a cutoff p±o ~ 2 GeV. This cutoff is also growing with energy, 
which makes it difficult to extrapolate safely to the high energies at LHC. 

The symmetry between the two ends of the parton chain implies that 
the LDC formalism also is applicable to hadron-hadron collisions'^. In the 
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Figure 4. Forward jet data from HI and ZEUS compared to the Cascade MC and LO 
and NLO dijet calculations^^. 

k_\_- factorization formalism the off-shell matrix element does not blow up 
when the exchanged transverse momentum k± 0, and we have therefore 
a dynamical cutoff for soft minijets. The cross section for a chain in pp 
collisions (which possibly may contain more than one hard subcollision) 
can thus be obtained from the fit to DIS data. An important point is here 
that the result is insensitive to the soft cutoff, Qq, in the evolution. DIS 
data can be fitted with different values for Qq, if the input distribution 
fo{x,Ql) is adjusted accordingly. If Qo is increased, the number of hard 
chains decreases, but at the same time the number of soft chains (for which 
all emissions have q±< Qo) increases, so that the total number of chains is 
approximately unchanged. 

There are two sources for multiple interactions: It is possible to have 
two hard scatterings in the same chain, and there may be more than one 
chain in a single event. The LDC model, when applied to pp collisions, 
can predict the correlations between hard scatterings within one chain, and 
also the average number of chains in a single event. The experimentally 
observed "pedestal effect" indicates that the hard subcoUisions arc highly 
correlated, so that central collisions have many minijets, while peripheral 
collisions have fewer minijets. In Pythia comparisons with data favour a 
distribution in the number of subcoUisions, which is very close to a geomet- 
ric distribution^^. 

Some preliminary results from the LDC model are shown in fig. 5. Here 
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a geometric distribution is assumed for the mimbcr of chains in one event. 
Fig. 5a shows the number of minijets in the "minimum azimuth region" 
60° < < 120° at = 1.8 TeV. The two LDC curves are obtained for 
soft cut-off values 0.99 and 1.3 GeV, showing the insensitivity to this cut- 
off. The two Pythia curves correspond to default parameter values, and 
parameters tuned to CDF data^^ . We note that the LDC result agrees very 
well with the tuned Pythia result. Fig. 5b shows corresponding results 
for LHC. Also here the two curves correspond to different cut-off values, 
and for comparison the result for 1.8 TeV is also indicated. We see that 
the activity increases by a little more than a factor of 2 between the two 
energies. 
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Figure 5. The average number of minijets in the "minimum azimuth region" for |?7| < 
2.5 vs. Ej_ for the hardest jet. Left: For ^ = 1.8 TeV. Right: For 14 TeV. 



The symmetry of the formalism implies that the chains join at one 
end at the same rate as they multiply at the other. The chain cross section 
grows like s^, and therefore the average chain multiplicity satisfies <nchain> 
oc .s^/(Ttot- Thus the results also may have implications for unitarization, 
saturation and diffraction. Work in these areas is in progress. 



5. Conclusions 

Our main conclusions can be summarized as follows: 

- Unintegrated parton densities are not observables, and their properties 
depend strongly on the definitions. 

- Different formalisms give similar results for J^{x, k'^,q'^ = 4:k^). This 
also implies that the predictions for observable quantities often are similar. 

- The role of the non-singular terms in the splitting functions is still a 
problem. 
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- Obscrvablcs without a large rapidity separation are often well de- 
scribed by higher order matrix elements plus DGLAP evolution. 

- There is a close relation between DIS and high energy pp collisions. 
The properties of the underlying event and minimum bias events can be 
predicted from DIS data. 
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